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^ . The frictional drag force between separated but coupled two-dimensional electron gases of dif- 

1 ^ ' ferent temperatures is studied using the non-equilibrium Green function method based on the sep- 
^— ( , aration of center-of-mass and relative dynamics of electrons. As the mechanisms of producing the 

frictional force we include the direct Coulomb interaction, the interaction mediated via virtual and 
' real TA and LA phonons, optic phonons, plasmons, and TA and LA phonon-electron collective 

modes. We found that, when the distance between the two electron gases is large, and at intermedi- 
ate temperature where plasmons and collective modes play the most important role in the frictional 
drag, the possibility of having a temperature difference between two subsystems modifies greatly 
the transresistivity. 
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I. INTRODUCTION 



Decades after the original proposal,0 and after the first successful measurement the frictional drag effect between 
separated but coupled electron subsystems is still attracting a lot of interest. In part this is due to the possibility of 
a direct view of interactions between charged carriers, which play a key role in central questions of condensed matter 
' physics likCrSjunerconductivity and quantum Hall effects. After the initial theories based on the directj-CnjjioBji^drag 
^ . mechanism,clfll3 electron-electron iiiteractions mediated by other channels, such as phonon exchangeaCl'EfLjIiJ'y and 
O ' electron-phonon collective modes,u have been explored as sources of frictional drag at low temperatures. However, 
, O, , recent experiments showed that further theoretical developments are still necessary before obtaining a quantitative 
view of this problenalij. At intermediate temperatures, the importpjce of plasmons has been predicted theoretically 
by Flensberg et aV0 and confirmed experimentally by Hill et aZ,ta but the discrepancies between theoretical and 
^ ' experimental results are significant. Giiven and TanatarEj made a detailed consideration of the optical phonon 
\l \ contribution to the plasmon drag, and obtained an improvement on the theoretical result for the case where optical 
2^ t phonons are of little importance on the drag force due to difficulty of exciting themtj at temperatures where plasmons 
domina te ., J he same time some studies concluded that the RPA should be substituted by more detailed many-body 
theoriestlnlSllj when estimating the frictional force. Besides those studies focusiiM|JMi-J|^j-a^ of frictiorial 



■ drag between 2D-2D electron-electron Fermi gas subsystems, intensive theoreticalt^lij|E2tEyE3cj and experimentalEll 
, works have been done on a 2D-2D electron-hole Fermi gas, and also in other case£a'E£lEZlG3. 

• Most works on Coulomb drag are made in the linear regime, in such a way that the temperature is assumed to be 

■ the same for the two subsystems, and equal ta-.that of the heat bath. It is well-known, however, that hot electron 
H ■ effects occur in semiconductor microstructuresEj, where electrons can be easily heated or cooled by very low electric 

\ voltage or photon excitations. In non-equilibrium situations, as it is the case in a Coulomb drag experiment, there 
is no reason to assume, from the very beginning, the two system at the same temperature of the heat b ath . . The 
question of thermalization of non-equilibrium systems is of interest also in several other fields in physicsEilEi. In 
^ this paper we study the frictional drag in the non-equilibrium regime resulting from the direct Coulomb interaction, 
;V ' the e-e interaction mediated by virtual and real TA and LA phonons, optic phonons, plasmons, and TA and LA 
phonon-electron collective modes. The frictional drag is described by a non-equilibrium Green's function method, 
which has been applied successfully in the development of quantum transport theory in semiconductors. We believe 
that allowing in the calculation different temperatures for the two systems, what is the central point in the present 
work, a new viewpoint is bcought to this developing area. 

In the Lei-Ting formalisnS used to describe non-linear transport in semiconductor structures, the balance equations 
are established by separating the center-of-mass from the relative motion. The Hamiltonian becomes H = Hq + He + 
Hi ^ He + He + Hr + Hi, where He corresponds to the motion of the center-of-mass. He corresponds to the 
time developing part due to external field, H^ is the Hamiltonian of relative electron system, and Hj represents the 
scattering by impurities and phonons. The last term is the one which causes the dissipation of momentum and energy. 

In order to obtain the equation of motion for the density matrix, we take as initial state an equilibrium state of a 
system free of the interactions Hj and He- The real system is obtained by adiabatically switching on Hj and He 
at the infinite past time, till the present non-equilibrium state. Then, the final state density matrix of the relative 
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electron system can be expanded in terms of the density matrix of the initial state. If we choose a proper initial 
state close to the final state of the relative electron system, a good result can be reached by keeping only the term 
in first order of perturbation. This approach has made the Lei- Ting balance equation an efficient tool to calculate 
the non-linear transport properties in semiconductors, since the numerical calculation involved is greatly simplified. 
However, this technique has not yet been employed in the electron-electron frictional drag phenomena, mostly because 
of the complexity of the many-body effect due to the inter-subsystem electron-electron interaction. We noticed that 
the Hi is important for configuring the ensembles of the system for this problem, so higher order terms must be 
considered. On the other hand, it is not convenient to include those higher order terms in the usual version of the 
formalism. In this paper, we propose a way to circumvent this difficulty in the case of the frictional drag problem. 

II. THEORY 

The Hamiltonian is written as h ~ Iiq + hj, with ensemble average density matrices g'^ and go corresponding to h 
and ho , respectively. The ensemble average of an explicitly time dependent Heisenberg operator At is: 

At g'^At >= Tr{goAt} ~ i ( dhTrig'^lA'l^ {t - ii),4°J}, (1) 

^ — oo 

where, 0^{t) — g^htQ^-iht ^^^^ 0^°{t) = e»''ot(9g-'''ot^ Instead of approximating g^ by go, as in the standard version 
of Lei-Ting's balance equations, we argue that, in cases where hi plays an important role in the ensemble averaging 
process, it is preferable to use A^ and h^ in place of A^° and /i'*". The validity of this approximation is the same as 
in the standard version, as we can see by considering the Hamiltonian ho ~ h — hj and expanding go in terms of g^. 
The average value of a physical property At becomes: 

At g'^At > = Tr{goAt} ~ i f diiTr{e''[A^(t - h), h'i^^] 

J — OO 

= Tr{BoAt}-if dt,Tr{g''[A'l{t),h%^{U)]} (2) 

Next, we consider a system of two coupled, but separated, quantum wells (no particle exchange between them being 
allowed), with carrier densities and electron temperatures ni and T^i in one well , and rt2 and Te2 in the other. As 
a general case, it is assumed that the two subsystems have drift velocities vi and V2, and applied electric fields Ei 
and E2, respectively. The temperature of the phonon reservoir is Tp. After separating the center-of-masses (located 
at positions Ri and R2) from the relative coordinates, the Hamiltonian becomes: 

H = eniRi-Ei + en2R2-E2 + i?ic + i^2C + H^i + H,,2 + i?ei2 + Hp,, + H^-ph + H^i, (3) 

where Hic = 5;^, = 2^' ^^i = Ek^k^^W + lJ2cl^Uc{ql)p\,plq,, H,2 - Ep '^p^^p + 

lEqs Uc{q2)pl^P-q.,, and iJei2 = Eqt^rf('?)PqP-qe"'^^^~^'^. The density operator in system 1 is p^ ^EkCk+qCk, 
and, for system 2, ~ Ek^iL+q^k- The final three terms, Hph, H^-ph, and H^i, are defined as usual.o Uc{q) and 
Ud{q) are the intra- and inter- well interactions. For the sake of simplifying the formalism, we start by considering 
the interaction only through the Coulomb electron-electron (e-e) potential. 

As in the original Lei-Ting balance equation formalism, we choose an initial equilibrium system composed of 
independent subsystems (electron gases in well 1 and 2, and phonons) in the Hamiltonian Hr — H^i + He2 + Hph, 
with the initial density matrix go = gei9e29p = pJ^p/'^pe^i2/Te2gHci/T^i ^ p^j, ^j^^ density matrix of the final relative 
electrons of the Hamiltonian Hr + Hj ^ Hr + Hei2 + H^-ph + Hei, we use g = e^f/'^" e^'=2/^''"+^'=i/^''i+^''i2/^/^^i^, 
which we assume to be close to the steady state of the relative electrons. The reason why we include i?ei2 into the 
density matrix is that the e-e interaction, which is very special and much stronger than other interactions such as 
electron-impurity and electron-phonon interaction, should give important contributions to the ensemble average for 
the final state of the relative electrons. The present choice accelerates the convergency. 

We consider the case where the electron temperature in well 2 is the same as that of the lattice, Tp = Te2 = T, 
and 7 — Tei/T. By introducing a new interaction representation corresponding to the virtual Hamiltonian H = 
'jHei + He2 + Hph + ^Hei2, where operators are expressed as 0{t) = e'^'Oe"*^*, the Liouville equation to the first 
order becomes: 
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g = go - I 



By approximating e 



iH{ti-t) ~iH(ti-t) 



-i(7-l)ffel(tl-t) 



the ensemble average becomes: 



At =< g^'At >= Tr{goAt} 



dtiOit ~ ti)Tr{g[A 



where Ay [t) = e'"^e~'^''-'^^"-^^ A^e'^''-'^''"-^^ e'''"K 

The time derivative operator of the total momentum for the electron gas in well 2 is: 

P2 - -i[P2, H] = -iVn.H = iqC/rf(g)pip2_qg-.q (Ri-R2) _ g^^Ea + P2.ep 



(4) 



(5) 



(6) 



In the case where the circuit of the second well is opened (without current), the phonon frictional force P2.ep and 
the impurity frictional force P2,ei disappear, together with the derivative of total moment of the second well. In the 
subsequent equations we make R2 = 0. The frictional force acting on the electron gas in the well 2 is balanced by the 
electric force: 



)]+iQ-Ri(ii)gi(l-7)(£fc+,-efc)(t-ti) 



/oo 
4 °° k,p Q 



(7) 



The above expression can be simplified if (i) the velocity fluctuations are neglected by making Ri(t) = v^r, where 
"Vd is the drift velocity of the electron gas in well 1, and (ii) we consider zero frequency drag, where the imaginary 
part of the frictional force is zero. Furthermore, the system satisfies the time invariance transformation in steady 
state with Q = 0. Then we have 



e7i2E2 = -/m ^ qr^(k, p, q,LJo) 

k,p,q 



(8) 



where ujq — qvd + (1 — 7)(£k+q — £k), and r^(k, p, q, w) is the Fourier transform of the retarded Green's function 

r«(k,p,q,T) = -*0(r)Tr {^[C/rf(<z)5[+q(r)5k(r)4_q(T)dp(T),C/d(-g)/5Lq(0)p2(0)]} (9) 



The later can be obtained by analytical continuation from its Matsubara counterpart: 



r^(k,p,q,L^) 



lim r'"' ik,p,q, iLUn) 



(10) 



After a series of calculation to expand r*^(k, p, q, ztJn) in the equilibrium representation without inter-subsystem 
interaction Hj, we obtain a diagram, where the renormalization is made with the polarization function of well 1 
and well 2, rather than inter-subsystem interaction. This diagram is similar to the result of the Memory Function 
method, but with different electron temperature in different subsystems. Following the suggestion of L. Zheng and A. 
H. MacDonaldjEl we adopt the renormalization of the Coulombic interaction between electron gases 1 and 2, as shown 
in Fig. 1. The drag resistivity then becomes: 



Rd - 



eniVd 



Tmin2e-^Vd 



E 



duj 



Ud{q) 



eT(q,t^) 



ni(q,c^,Ti)n2(g,cc> + Cc>o,T) 



(11) 



The drift velocity of electrons in well 1 is assumed along the x direction, q is the carrier momentum, uj is the 
carrier energy, n{x) is the Boson-Einstein distribution function, n^(q,w,Ti) — Il}(q,a;,Ti) ~\- iti\((\,ijj,Ti), and 
n^(q, a;,T) = II^(q, WjT) -|- illKq, Ci^, T), are the zero-order polarization functions of carriers in well 1 and well 2, 
respectively. The dielectric function of the electron system is 
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eT(q, t^) = [1 - f/c(g)ni(q, c^, ri)][l - Uciq)fi'iq, uj, T^)] - Uj{q)U\ci, to, Ti)n2(q, lu, T2) 



(12) 



These expressions were obtained by assuming that only the Coulomb interaction is important to the electronic 
coupling. It is easy to extend this result to the case in which other mechanisms of interaction between electrons 
become important, just by replacing the inter- and intra-well e-e Coulomb potential by the sum of the corresponding 
interaction potentials due to different mechanisms. For instance, if phonon mediated electron-electron interaction is 
to be taken into account, the following substitution should be made: 



Uc{q) 
Udiq) 



C/i2(q,w) 



Uciq) 
-- Ud{q) 



A 

E 

A 



(13) 



Here D^^ = / ^|M^(q, c.)| V.(g.)/,(g.)[- 



is the interaction potential between electrons 



in wells i and j mediated by phonons of the branch A, as discussed in Ref. 0. M'^(q, uj) is the bulk matrix of electron- 
phonon interaction, and fi{qz) is the form factor for the electron-phonon interaction in well i. The mean free time 
of phonons r^^ equals the mean free path divided by phonon velocity v\. For small drift velocities, the drag 
resistivity becomes: 



Rd — 



1 



n/2 



j fii(q,c^,r, 



d9 cos 



q^dq 



duj 



eT(q, w) 



n2(q,^,r) 



d 



T2Sm]i'{Lu/2T2) 



2 [coth(tj/2ri) - coth(w/2r2)] — n^(q, w,r2) 



(14) 



The real (virtual) phonon contribution comes from the imaginary (real) part of D,^- and the contribution of plasmons 
and electron-phonon collective modes come from the zeroes of the real part of eT(q, w). The well-know expression for 
the linear drag resistivity is recovered if we neglect the anisotropy of the interaction matrix, and assume the same 
temperature Ti = T2 = T for the electron gases in both wells. Then, the drag resistivity of the two-temperature 
model becomes the generally accepted linear result: 



Rd 



1 



87r2e-^nin2T Jq 



dqq^ 



duj 



Ui2(q, w) 



smh'^{uj/2T) '' 



(15) 



III. RESULTS AND CONCLUSION 



The numerical calculation is performed in a model in which the electron system occupies two GaAs quantum wells 
separated by infinite barriers. This model may underestimate the electron effective mass and the correlation between 
two wells, but these effects are negligible for the cases discussed in this paper. In addition, it is well known that 
the RPA is expected to give poor results for low density electron gases, and several improvements beyond RPA have 
been employed in the literature, as already discussed in Chap. |. However, RPA is used in this paper for estimating 
the many-body effects, since the mechanism leading to the unexpected large measured drag resistivity has not yet 
been clarified, and also because our aim here is not focused on the improvements or more detailed estimative or 
corrections beyond RPA. We believe that local field corrections to our result give similar results as discussed in other 
places. It is worthwhile to mention that effects of the optic phonon directly causing, or mediating, the frictional 
drag were considered in Refs. 15 and 0. Here the optic phonon cpBtribution is taken into account just by using the 



static dielectric constant Ks as the background dielectric function.Eij An argument in favor of this assumption is that 
the plasmon energy is much smaller than those of optic phonons, and at the plasma frequency the optic phonon 
contribution to dielectric function is practically frequency independent. 

The contributions of the longitudinal acoustic (LA) and the transversal acoustic (TA) phonons via deformation 
potential and piezoelectric interaction, on the other hand, are fully included. The mean free times, TiOp/j and Tropfi, 
are used as fitting parameters for both the LA and the TA branches. The parameters used for GaAs are Kg = 12.9, 
di = 5.31g/cm^, E = UeV, m* = 0.067, v^i = 5.29 x lO^m/s, cm = 1.41 x lO^V/m, and Vst = 2.48 x 10^ for the 
dielectric constant, density, deformation potential, electron effective mass, longitudinal sound velocity, piezoelectric 
constant and transversal sound velocity, respectively. The numerical calculation is performed for a structure of two 
200A wide quantum wells, with a 300A barrier between them. 
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In Fig. 2 (a), the drag resistivity Rd is plotted as a function of the equiUbrium temperature in (Ti = T2) for 
a matching electron densities of n = 1.37 x lO^^cm"^ in both wells. The calculation is performed in a range of 
values of r^^ and r^^ suitable for analyzing the relative importance of the electron-phonon collective mode for the 
LA and TA phonons. The solid lines represent results of different values of the phonon mean free paths (the same 
for LA and TA phonons), namely Iph — 10^^, 1, lO^cm. In the two upper curves (lcm,10^cm) the electron-phonon 
collective modes dominate the contribution. In the lower curve (IXt^cm), however the main contribution comes from 
the virtual and real phonon terms. Similarly to previous studieaS'til, our theoretical results fit qualitatively well the 
experimental data from Ref. |l^ (filled circles), but several quantitative inconsistencies still remain. For instance, at 
low temperatures the calculated peak of the acoustic phonon is displaced to the high temperature side. If we allow 
Ila and It a to be different and diminish the value of Ila to reduce the contribution of LA phonons, the peak moves 
to the low temperature side. In dotted line, we show the drag resistivity including only the contributions of the TA 
phonons {It A = lO^cm), together with the Coulomb interaction. Comparing with the experimental result, this seems 
to indicate that TA phonons have a longer mean free path, then LA phonon, which can result in the dominance of TA 
phonons in the frictional drag force as mentioned in RefJl^. However, even though the position of the calculated peak 
can coincide with the experimentally observed by fitting Ila and It A, the value of the maximum cannot reach the 
experimental result in this low carrier density case. At intermediate temperatures (T ^ 0.2Tf), the plasmon begin to 
dominatte the contribution to frictional drag force, which reaches its maximum at about (T ~ O.bTp), as calculated 



in Ref.ta. Our conclusion is t 



it the LA phonons enhance the plasmon peak at intermediate temperatures, in the 



IS t lat 

same way as the optic phonon,E£l, while the TA phonons have little influence at this temperature range. However, our 
calculation gives a smaller frictional drag force, and a higher peak temperature than what is observed experimentally. 
In Fig. 2(b), we show the dependence of the drag resistivity on the temperature, as in Fig. 2(a), with a different carrier 
density, n = 2.66 x lO^^cm"^, and several TA and LA mean free paths, namely, Iph = 10~^, 10~^, 1cm. Differently 
from the previous low density regime, we find a very short phonon mean free path, Iph = 0.01cm, when fitting the 
experimental result. Since samples differing just by the carrier densities are expected to have the same order of the 
phonon mean free path, this points to a disagreement between theoretical and experimental results for the density 
dependence of the frictional drag force. If the LA phonon is excluded from the calculation, as denoted by the dotted 
curve, even a very long mean free path for TA phonon, Ita — lOOcm, results in a very small frictional drag force. 

Next, we relax the condition Ti ^ T2, allowing for the driving well to have a different temperature from that of the 
drag well. In Fig. 3(a) and (b) we calculate the dependence of the drag resistivity on the normalized temperature 
divided by the Fermi temperature of the drag well, T2/TP. Electrons in the drag well are assumed in thermal 
equilibrium with the phonons. Curves are plotted for different ratios T1/T2 ■ It is found that a small change of that 
temperature ratio affects little the phonon drag peak at low temperatures, but modifies greatly the plasmon peak at 
intermediate temperatures. When the driving electrons are cooler than the drag electrons, the frictional drag force 
increases greatly. In the reverse case the plasmon peak diminishes, or even disappears. If we take a detailed look at 
the frictional drag, given by Eq.(I^, we find that the changes introduced by non-equilibrium come mainly from the 
second term. It reflects the modification of the correlation function by the drift energy of electrons in the driving 
well, and it is of second order in the equilibrium. 

In Fig. 4(a) we calculate the non-equilibrium drag resistivity divided by that obtained in equilibrium, as a function 
of the temperature ratio T1/T2, at various T2, in the case of matching densities ni — n2 — 2.66 x lO^^cm"^. At low 
temperatures (S-fiT), a small change on the temperature of the driving electron gas makes little effect on the value 
RdT2/ R^Ti, as shown by the almost flat solid curve at RdT2/R^Ti = 1. At higher temperatures, however, when the 
interaction with plasmons becomes important, a cooler (hotter) driving electron gas enhances (reduces) greatly the 
drag resistivity. The derivatives of curves in Fig. 4(a) at the point Ti = T2 are plotted in Fig. 4(b) as functions of 
the temperature of the drag well. Results are also shown for the matching electron density n = 1.37 x lO^^cm"^. The 
sensitivity of the plasmon contribution to the temperature difference ST = Ti — T2 is a direct result of the non-linearity 
of the correlation function over frequency, conflrming that a correct approach to the Coulomb drag effect must contain 
the non-equilibrium assertive of the present work. 

In order to make the above point clearer, we obtain, from Eq. (|l^), the expression for the resistivity difference SRd 
resulting from a temperature difference near T2 — T. 

fdqq^ rd.|W,.(q,.)|^5M^i;;^^ (16) 
T 87r^e^nin2T Jq Jq smh (w/2r) 



where Q(q,u;) — — - -^^Iliiq,^) ■ If H2(q,ijj) is a linear function of uj, then Q(q,u;) = 1 and SRd = -TpRd, 
i.e. Rd/{TiT2) is a constant near Ti = T2. Otherwise, Rd/{TiT2) will be sensitive to the temperature difference. The 
contribution from acoustic phonons comes mainly from q ~ ikp, where Il2('?, w) is a linear function of w in the range 
of the corresponding phonon energy (with the parameters used in this paper) . This leads to the fact that the acoustic 
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peak in Fig. 3 does not show any non-equilibrium trace. On the other hand, the plasmon energy is higher than the 
energy of the single particle excitation, and the corresponding non-linearity of ^) invokes strong non-equilibrium 
character for drag resistivity Rd/{TiT2) due to plasmons. 

Fig. 5 shows the drag resistivity Rd/{TiT2) in equilibrium cases (solid curves) as a function of the ratio ni/n2 of 
the electron densities at different drag temperatures, T2 = 0.04, 0.09, 0.16, 0.24, and 0.44Tf in comparison with the 
experimental results. The mean free paths of the LA and TA phonons is assumed to be Iph = O.Olcm. We found that 
the acoustic phonon enhances the drag resistivity at large ni/n2. It can be expected that the inclusion of local field 
effects will provide good fittings between theoretical and experimental results,EJ except at temperature T2 = OAATp- 
The results with bare Coulomb interaction is shown in equilibrium cases by dotted curves, and in non-equilibrium 
cases with Ti = O.9T2 by dashed curves. 

In conclusion, we have studied the momentum transfer rate between two nearby separated electron gases, which 
are coupled via the Coulomb interaction and phonon mediated interaction. We focused on the non-equilibrium 
configuration, and a theory was developed to describe the frictional drag force felt by one electron gas as a result 
of the relative drift of the other, taking into account the possibility of difference on the subsystems temperatures. 
It was found that a cooler (hotter) driving electron gas greatly enhances (decreases) the frictional force caused by 
plasmons. This behavior results from the fact the plasmons locate in a region of the lo — k space where the correlation 
function shows a nonlinear frequency dependence, and they modifies little the acoustic phonons force because the 
latter are important mainly in the region of the energy-momentum space where the correlation function is linear with 
the frequency. 
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FIG. 1. Feymann diagram for the drag resistivity 

FIG. 2. Drag resistivity Rd/T^ as functions of temperature T /Tf at electron density (a) n = 1.37 x lO^^cm"^ (Tf — 57K), 
and (b) n = 2.66 x lO^^cm"^ (Tf ~ IWK). Results including the direct Coulomb interaction plus the LA and TA phonons 
with the same mean free paths Ip^ are plotted by solid lines. Experimental results according to Ref. |l^ are shown as filled 
circles. The dashed line shows the results when LA phonons are excluded in the calculation to better identify the peak position 
of TA phonon. The phonon mean free paths are denoted by numbers besides the curves in cm. 

FIG. 3. Drag resistivity as functions of drag temperature T2 at different driving/drag temperature ratio T1/T2, denoted by 
numbers besides the curves. The electron densities and phonon mean freepaths are (a) n = 1.37 x lO^^cm" , Iph = 100cm and 
(b) n = 2.66 x lO^^cm"^, Iph = 0.01cm. Experimental results from Ref. |lj are swhon as filled circles. 



FIG. 4. (a) The normalized drag resistivity RdT2/R^Ti as functions of driving/drag temperature ratio T\/T2 at fixed drag 
temperature T2 numbered besides curves, and (b) its derivative at T\/T2 — 1. Only the direct Coulomb interaction is considered. 



FIG. 5. Drag resistivity Rd/{T\T2) is shown as function of driving/drag electron density n\/n2 at different normalized drag 
temperature T2/TF = 0.04, 0.09, 0.16, 0.24, 0.44 and driving electron density m = 2.66 x 10^^ cm~^. Experimental results from 
Ref.lij is shown by filled circles. Solid lines represent results including contributions from all kinds of mechanisms in equilibrium, 
with TA and LA phonon mean free paths Iph = 0.01cm. Dotted lines represents direct Coulomb drag forces (including plasmons) 
in equilibrium. Dashed lines show the Coulomb drag forces in non-equilibrium cases with T\/T2 = 0.9 
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